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Abstract 



The MEDLEY setup based at The Svedberg Laboratory (TSL), Uppsala, Sweden has previously been used to 
measure double-differential cross sections for elastic nd scattering, as well as light ion production reactions for 
various nuclei in the interaction with neutrons around 95 MeV. When moved to the new beam line, the first 
experimental campaign was on light-ion production from Ca at 94 MeV in February 2005. 

These data sets have been analyzed for proton production in forward and backward angles. The Ai? — Ai5 — 
E technique have been used to identify protons, and a cutoff as low as 2.5 MeV is achieved. Suppression of 
events induced by neutrons in the low-energy tail of the neutron field is achieved by time-of-flight techniques. 
The data are normalized relative to elastic np scattering measured in the 20-degree telescope. Results from an 
accepted neutron spectrum are presented and some methods to correct for events from low energy neutrons are 
presented and evaluated. The data are compared with calculations using the nuclear code TALYS. It was found 
that TALYS systematically overestimates the compound part, and underestimates the pre-equilibrium part of the 
cross section. 
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1 Introduction 



In the last years, fast neutrons have gathered more and more interest worldwide. Lots of applications involving 
fast neutrons in some way are predicted to be available within 50 years, for example one can mention fusion 
power reactors and transmutation of spent nuclear fuel. Even today, fast neutrons play an important role in the 
field of electronics, where it can cause so called single event effects, doing damage to the hardware and software. 
Fast neutrons are also used for hadronic radiotherapy of cancer in several places in the world, for exampel at 
iTemba LABS in South Africa. So understanding of fast neutron interactions is a large and important part of 
applied nuclear physics. 

In this work the interaction between 94 MeV neutrons with calcium will be examined, in particular the 
productions of protons at forward and backward angles. A first rough estimation of the angular distributions 
will also be carried out. In this introduction to the work, some general properties of the neutron, and experiments 
with the neutron, will be discussed, as well as a short introduction to the nucleus under study. This is followed 
by a few sections about what is actually measured and some possible applications and motivations for this work 
in the form of a biological introduction. 

1.1 The neutron 

The second lightest (known) baryon in the universe is the uncharged neutron, discovered in 1932 by James 
Chadwick 1 1], when he was examining mysterious radiation emitted from beryllium, boron and lithium that had 
been discovered a few years earlier by Bothe and Becker |2 1. Since the neutron is the second lightest, the slight 
mass difference between it and the lightest, the proton, makes the neutron unstable when it is in an unbound 
state and it decays according to figure [TT| with a lifetime of 887 seconds |3J. 



e 
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Figure 1.1: fi-decay of a free neutron into a proton, an electron and an antineutrino. 

The neutron interacts with its surroundings through all four of the fundamental forces: the strong nuclear 
forc^ the weak nuclear forc^ the electromagnetic forc^ and the gravitational forc^ In experiments with 
neutrons there exist a few difficulties that do occur due to its special properties. Acceleration of neutrons is 
impossible today since accelerators use the charge of the particles. Instead, one has to rely on a neutron beam 
as a secondary beam. This means that one uses a primary beam of charged particles, for example protons, and 
let them hit some target that can induce a neutron production nuclear reaction. Similarly the focusing of particle 
beams is done by large magnets, while the beam shaping for a beam of uncharged particles must be done with 
collimators. Dealing with neutrons also gives rise to background and radiation protection issues. This since 
high-energy neutrons has a lower probability than charged particles for interaction and higher possibiUty for 
passing through the shielding material. 

'More or less all nuclear reactio ns. 
^In the /3-decay showed in figure [LT| 

'Even though the neutron is electrically neutral it still holds a magnetic dipole moment of about -1.91 nuclear magnetons, /in- This is 
due to the fact that a baryon has a magnetic moment equal to the sums of the quarks magnetic moments, which in the case of the neutron 
gives Un = I A^d — \t^a = — % — — ^ ■ The electrical dipole moment is measured to be less than 2.9 ■ 10~^® e-cm |'4']. 

''For example in neutron stars. 
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1.2 Calcium 



The nucleus under study in this work builds up a quite hard, silvery white, alkaline earth metal. The name 
calcium origins from the Latin word calcis that means lime, referring to the group of minerals dominated by 
calcium-carbonates, -oxides and -hydroxides often used as building material, for example in Umestone, concrete 
and mortar between bricks in masonry fSl . 

In the human body, calcium is the most common mineral where almost all of it is located in the skeleton in 
the form of hydroxide apatite crystals. The skeleton of a normal adult contains about 1 kg of calcium. Some of 
it can also be found in the soft tissue in the form of Ca^+ ions f6l. There are six stable isotopes of calcium where 
the doubly magicP] isotope "^°Ca is the absolutely most common one, with a natural abundance of 96.941 % ||3|. 



1.3 Cross-section measurements 

One of the most fundamental needs in nuclear physics, as well as in all other physics and even all other sciences, 
is theoretical models to explain nature. Of course, another fundamental need in the field of science is also the 
to go the other way, by obtaining experimental data to verify theoretical models. 

In nuclear and particle physics the cross section is in a way the connection between everyday life and the 
often complicated and abstract theoretical models. It is a quantity that both can be measured in experiments, as 
well as calculated from simple or complex nuclear models. The unit is a non-SI, although accepted for use with 
the SI, areal unit named barrj^that can be visualized as the effective area a projectile sees when approaching a 
target. This makes the interpretation of a cross section as a reaction probability a little more intuitive; the larger 
the effective area of the target - the larger the probability of hitting it. It is important to keep in mind though 
that the reaction cross section is not the same as the geometric cross section. 



1.4 Ionizing radiation and humans 

The human body is constantly exposed to radioactivity in terms of ionizing radiation. The sources are partly due 
to human activities like rest products from atmospheric nuclear weapons tests, nuclear power accidents, dental 
x-rays and many other things. But there are also various natural sources of radioactivity. For example radon gas 
from earths crust and cosmic radiation. Even man itself contains various radioactive isotopes and is a source for 
radiation, as illustrated in figure [TT2| 

The damage that ionizing radiation does on human tissue is mainly on the DNA and a majority of the DNA 
damage is indirect by ionizing the water in the cell and creating very reactive free radicals. Following |9|, there 
are three possible damages that can occur, namely base pair damage, single strand break and double strand 
break. Base pair damage and single strand break are the most common of these three, but causes generally no 
problems for the cell and are fixed fast and easy. Double strand break on the other hand can create some more 
severe problems. 

When the DNA, in the form of chromosomes, gets exposed to a double strand break and the cell tries to repair 
it there are three possible outcomes. The first is that the reparation succeeds and everything is back to normal, but 
the reparation can also go wrong in two different ways. One way is that we get an abnormal chromosome, like a 
ring chromosome where the end of each arm has been lost and repaired into a ring formation. Another abnormal 
chromosome, that will appear in about 60 % of chromosome aberrations, is the di-centric chromosome where 
the chromosome has been repaired into one with two centromeres, the region where the two the two identical 
parts of the chromosome touch during cell division, instead of just one. Both these examples of abnormal 
chromosomes are, although fatal to the cell, quite harmless to humans since the damaged cell never can survive 



mitosis, cell division. This is illustrated in figure 1.3 

The third outcome, also illustrated in figure |1.3| of a double strand break is the worst one. Here the chro- 
mosomes undergo stable rearrangements instead of getting abnormal shapes. The stable rearrangements could 



See section 



2.1 



^Rumors claim that Baker and HoUoway 17 1 coined the term bam in December 1942. Trying to come up with a name for the unit 10 

m^ 

in the experiment was as high as the probability of hitting a barn when firing a gun. 
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Figure 1.2: In some sense, everyone of us is a Radioactive Man j&l/. 



consist of translocation of two of the arms or that a small portion of an arm is removed from the chromosome. 
This mutation does not lead to mitotic cell death. Instead the damage will multiply at the same rate as the cell, 
something that very well can be a part of inducing a cancer tumor. This is the reason why cancer-inducing 
effects are often refereed to as stochastic effects. The probability of the statistically rare event of mutation 
increases with increasing exposure, but the effect of an injury is independent of the received dose. 
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Figure 1.3: Possible reparation errors when fixing a double strand break and the resulting DNA-damage flO^. 



1.5 Cancer therapy 

In principle there are three ways to treat cancer. One can use surgery, chemotherapy or radiation therapy. The 
basics behind radiation therapy of cancer are to damage the cancer cells DNA strands and in that way prevent 
the cell to divide. In this case one actually wants a double strand break that ends up in an abnormal chromosome 



and prevents the mitosis, as mentioned in section 1.4 The most common type of radiation therapy is using 7- 
rays from ^"Co, x-rays or electron beams. However some tumors, called radioresistant tumors, seem to respond 
quite poorly to this treatment and therefore other methods are needed. 

One way is to use hadrons such as protons, neutrons, helium ions or pions instead of photons or electrons. 



Since hadrons have a larger linear-energy-transfer ( LET 1, that is the stopping power in the medium, than elec- 
trons and photons, they has a higher chance to do biological damage to the cell. Actually there are three things 
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that make a high LET particle more effective than a low LET particle. Firstly it is harder for the cells to repair 
the larger damage, secondly high |LET] particles is more efficient in damaging the oxygen-poor tissue in a tumor. 



And finally high LET radiation is not as sensitive to whether or not the cell is currently dividing as low LET 
radiation f9"TTl. Another property of charged hadrons is the interaction depth. While photons and neutrons 
deposit most of their energy in the surface, charged hadrons deposit most of their energy in a well-defined depth 
in the patient. This known as the Bragg peak. 

Using fast neutrons for cancer therapy is a little more complex method than using photons, electrons or 
hadrons, since neutrons do not damage the tissue via direct interactions. Instead, nuclear reactions with neutrons 
damage the tissue indirectly by production of various types of light ions. This means that in order to make 
accurate dose calculations, typically via Monte Carlo based radiation transport codes, it is fundamental that the 
double-differential cross sections for light ion production are well known lil2J . 



1.6 TALYS 

TALYS is a joint project between NRG Petten in the Netherlands and CEA Bruyeres-le-Chetel in France, and 
is a nuclear reactions modelling program working with nuclear reactions of energies 1 keV - 200 MeV. It is 
written completely from scratch in Fortran??, with the only exception of the implemented coupled-channels 
code ECIS-9? ifTSll . 

The two main purposes of TALYS are to be used as a physics tool during experimental analysis, and as a 
data tool to generate nuclear data for use in various applications. The program is also very flexible, since one 
can get a reasonable cross section by a minimal input like the four lines: 

projectile n 
element Ca 
mass 4 

energy 94 

But it also reacts to more than 150 keywords, making it possible to adjust the nuclear models and the output 

m. 
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2 Theory 



As mentioned in section 1.6 the theoretical predictions are made using a nuclear reactions code called TALYS. 
This is an outline of how TALYS does the calculations in this particular case, after a brief introduction to the 
basics of nuclear reactions. 

One should remember that this is only a crude simplification of the calculation process valid only for inci- 
dent high-energy neutrons on a spherical nucleus and emission of protons. When dealing with more complex 
particles like deuterons, tritons, helium isotopes and so on effects like continuum stripping, pick-up and knock- 
out also needs to be taken into account. Deformed nuclei also introduce a number of corrections not needed 



here. An example of the different reactions contribution to the total cross section can be found in figure A. 1 



2.1 The shell model and magic nuclei 

In nuclear physics there is not a single simple model to describe the atomic nucleus completely. Instead several 
different models that fit different kinds of parameters or properties for the nuclei are used. For example there 
are the liquid drop model and the shell model. 

The shell model for nuclei is quite analogous to the shell model for atoms. In the atomic shell model the 
electrons populate shells outside the nucleus and when a shell is filled they start filling the next. In a nucleus 
the principle is the same, only here it is protons and neutrons that occupy shells and not electrons. Since the 
neutron and the proton have individual shells two layers of shells are filled in parallel [IS]. 

Since the shell model deals with particles populating different energy levels there are some regions of extra 
interest and these are the nuclei with filled, or close to filled, shells. The nuclei with certain number of nucleons, 
corresponding to filled shells in the shell model, are found to be more tightly bound than others. These numbers 
of nucleons are commonly referred to as magic numbers and are: 2, 8, 20, 28, 50, 82 and 126. Nuclei containing 
filled shells are consequently called magic nuclei and nuclei with both protons and neutrons in filled shells are 
called double magic nuclei. Examples of double magic nuclei are ^^O, ^°Ca, **Ca and ^°®Pb. 

There are theories that also higher numbers of neutrons and of protons are magic numbers which would give 
rise to stable exotic nuclei such as ^?^4Uuq, f^pUbn and fjeUbh in an island of stability among the otherwise 
highly unstable elements at the end of the periodic table [161 . 

2.2 Fundamental nuclear reactions 

One important field of research to understand nuclear physics is the study of nuclear reactions. The general 
nuclear reaction is when an incoming particle, a, hits a target, X, causing the emission of a particle, b, and a 
recoil nucleus, Y. This can be written on the same form as a chemical reaction 

a + X^Y + b (2.1) 

or in a more compact form X{a, b)Y where the Y often is dropped. An example|^of a reaction could be 

p + ^Li ^ *He + a (2.2) 



or correspondingly '^Li(p, a). One can now divide (2.1 1 into two subgroups, elastic and inelastic nuclear reac 



tions. An elastic reaction is when a — b and X — y in (2.1 1 and no energy is lost due to excitation of any of 
the components, while in an inelastic process the target gets excited and/or a ^ b. 

When studying high-energy reactions, another way to classify nuclear reactions is in terms of the time scales 
at which they occur, or equivalently in terms of the number of collisions within the nucleus. The fastest reac- 
tions, with the fewest collisions, are called direct reactions and occur on time scales of about 10^^^ s. Here the 
incoming projectile primarily interacts with the surface of the target. At the opposite end there is the compound 
reaction at time scales of about 10~^^ — 10~^* s 1, 15 1 , and in between these types pre-equilibrium reactions may 



occur. The spectrum of the different reactions is illustrated in figure 2. 1 where also the corresponding shape of 
the cross section is sketched. 



'This example is actually quite historically interesting since it was the first accelerator induced nuclear reaction fl51 
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Figure 2.1: The different types of reactions, adapted from HMV : Compound reactions are indicated by C, pre-equilibrium 
reactions by P and direct reactions by D. The x-axis represent reaction time and particle energy and the y-axis the cross 
section. 



2.3 The optical model 

Before diving into details in the different nuclear reaction it is suitable to have a closer look at the basis of nuclear 
reaction calculations, something that is called the optical model potential ( |OMP| i. An analogy to this nuclear 
model is a murky glass sphere with some incoming light. Some of the light will be refracted, or elastically 
scattered, through the sphere, while other rays will be partially absorbed. In the nuclear model, this effect is 
mathematically described by defining a scattering potential as 

U{r,E)^V{r,E)+iW{r,E) (2.3) 

that can be inserted into the Schrodinger equation. Following the parametrisation in lfT4l[T7l . these potentials 
can in turn be described as a sum of potentials caused by different properties of the nucleus as 

V(r, E) = -Vy(r, E) + Vso{r, E)-\.-a + Vc(r) (2.4a) 

W{r, E) = -Wy(r, E) ~ Woir, E) + Wso{r, E)-l.-a (2.4b) 

where the indices denote the volume part, V, the spin-orbit coupling part, SO, the surface derivative part, D, 
and the Coulomb part, C. Finally, all of these can be further broken down into an energy dependent part and a 
radial dependent form function part as follows 

Vy(r, E) = Vv{E)f{r, Ry, ay) (2.5a) 

2 . . 

Vso{r,E) = Vso{E) [ — ] -- f{r,Rso,aso) (2.5b) 

- r ar 

for r < Rn 

- (2.5c) 
for r > Rc 

Wv{E)f{r,Rv,av) (2.5d) 

Woir^E) = ^AaDWDiE)^f{r,RD,an) (2.5e) 

ar 

Wso{r,E) = Wso{E)( --^f{r,Rso,aso)- (2.5f) 

The parameters Ri = ViA^I^ and are the radius and diffuseness parameters for the nucleus and the form 
function often is assumed to be of Woods-Saxon type 

f{r,Ri,ai)^ ^T^H-- (2-6) 
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These parameters are listed by Koning and Delaroche fTT^ for a vast amount of nuclei and for **°Ca and 

The change in the potentials with energy of the incoming 



incoming neutrons they are given by table 2.1 



neutron and distance from the center of the nucleus is plotted in figure [2l2} 
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Table 2.1: Numerical values for the OMP parameters for Ca with an incident neutron. 




20 40 60 80 100 120 140 160 180 200 
Incoming neutron energy [MeV] 




Figure 2.2: Optical potentials plotted for '^'^Ca and incoming neutrons. The left panel shows Vv{E), Vso{E), Wv{E), 
Wd (E) and Wso (E) for incident neutron energies between and 200 MeV while the right panel shows V(r) (solid line) 
and W{r) (dashed line) for an energy of 94.5 MeV. Parametrized via HT^ . 



2.4 Direct reactions 



To calculate the cross sections from the direct formation processes TALYS has several models incorporated, 
Uke the weak coupling model or giant resonances. For near- spherical nuclides TALYS calculates the cross 



sections for the direct formation processes with the distorted wave Born approximation (DWBA i. This basically 
means that one uses first order perturbation theory on the regular Bom approximation, that the incoming and 
outgoing wave functions are just plane waves. This will give a result where the simple plane waves of the 



Born approximation is disturbed by the scattering potential of Eq. (2.3i. This first order model for smal l 
deformations is valid since "'"Ca is a double magic nucleus, and hence near spherical. For a reaction like (2.1 
the cross section is obtained from the relation 



do-" ^ EgEb Pb 1 |/f|/T r\|^ 



(2.7) 



where Ei is the total energy, pi the momentum, Si the spin of the particle, i, and nii the mass. The initial and 



final states (/| and \i) are given by the DWBA and Hint is the interaction potential lITSI . 
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2.5 Pre-equilibrium reactions 

After the direct reaction and before the nucleus reaches equihbrium some pre-equihbrium processes takes place. 
In TALYS the pre-equilibrium stage is calculated using a model called the two-component exciton model, 
following lfT4l \T9l . In analogy to solid state physics one can define particle-hole pairs, also called pairs of 
excitons, above the Fermi surface^as 



^The Fermi surface is the highest energy level when all the nucleons in a nucleus, or electrons in an atom, are at their lowest lying states 
- the Fermi sea. 
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(2.8a) 



(2.8b) 



h = h^r + h,. 



(2.8c) 
(2.8d) 
(2.8e) 



where p is the number of particles, h the number of holes and n the number of excitons while tt corresponds to 
protons and v to neutrons. When the incoming projectile interacts with the nucleus it can excite a particle-hole 
pair On its path through the nucleus it can either excite new particle-hole pairs, annihilate existing ones or 



change configuration of the existing pairs as illustrated in figure 2.3 as well as release some excitation energy 
by emitting a particle. The more interactions the particle undergoes on its way through the nucleus, the more it 




,0,3,3 



4,3,0,0 



Figure 2.3: Different possible reaction paths in the two component exciton model for an incoming proton [ 19]. The numbers 
in the boxes correspond to p^, /i^, and of some population P (p^, /i^). The thickness of the arrows is 
proportional to the transition rate, and as the arrows pointing backwards are really small compared to the others, they are 
neglected by TALYS. This is the so called never-come-back approximation. 

looses its memory of incoming direction and energy until it reaches some equilibrium state and the compound 
model takes over. The cross section for this pre-equilibrium state is described by 



do-E ^ CF 



E 



E 



Wk {Pn:K Ek) T (p^, /l^,p^, /ij,) P {pT,,h^,p,y, 



(2.9) 



where a'"^ is the total compound formation cross section, Wk is the emission rate, r the lifetime of the particle- 
hole states and P the states that survived and are left from the last process. 

The to tal compound-formation cross section is just the reaction cross section, t hat can b e calculated from 
the OMP with the sum of the cross sections for direct reactions calculated from the DWBA subtracted, a'^^ ~ 
^reac _ g-duect -pj^g emission rate depends on known parameters, the inverse reaction cross section calculated 



from the OMP and the particle-hole state density. Both the lifetime and the survived population are calculated 



from the internal transition rates for particle-hole pair creation, conversion and annihilation 1741 . 



2.6 Compound formation 

For the compound formation one can separate two different processes. At low energies it is a binary reaction 
where the projectile is captured by the target, followed by emission of another particle. But at high incident 
energies the nucleus will still be left in a highly excited state after the binary reaction and other process take 
over. Although the two processes at first sight seem identical, there are two major differences. One of them 
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being that the first process is non-isotropic, and the other the presence of a width fluctuation factor in the first 
process that correlates the incoming and outgoing waves [13 \. An analogy is often made with the evaporation 
of a hot fluid, that works in basically the same way, where particles evaporate from the nucleus until it is cooled 
down. 

Also for compound reactions TALYS includes several models for different situations. At high energies, as in 
this case, the default model used is multiple pre-equilibrium emission within the exciton model. This is an exten- 
sion of the pre-equilibrium model and includes the pre-equilibrium population {Z, N, pTr,h.^,p^,h^, (z) ) 
that needs to be included in (2.9 1 and summed over all possible combinations of parameters. The populations 
that are not used for a new V'^'"^ are used to feed another multiple emission model called multiple Hauser- 
Feshbach decay lfT4ll . 



2.7 Angular distribution 

Continuum angular distributions can be achieved using the systematics by Kalbach EOl . The systematics is 
based on that the angular distribution, to a first approximation, is independent of the energy of the projectile, 
as well as the type of the target, projectile and ejectile. Instead it depends on the emitted particles energy 
and the fraction of multi-step direct and multi-step compound emissions for the given energy li20J . In this 
phenomenological approach the angular distribution is described as 



d^cr 1 da 



AnAE Att dE sinh(a) 



(cosh(acos0) + /MSDsinh(acos0)) (2.10) 



with a = a{E) as a free parameter to be fitted and /msd as the fraction of emissions that comes from multi-step 
direct processes. This parameter, as well as the energy-differential cross section, 4?;, is assumed to be known. 



This result is physically motivated by Chadwick and Oblozinsky 1.211 . The form in Eq. (2.10 1 is the form that 
TALYS uses, but for fitting purposes the generalized form 



dV 



= 6e°™"^ (2.11) 



AVLdE 

is used. Here a = a{E) is still a free parameter and another free parameter, b ~ b{E), is introduced. 
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3 Experimental setup 



The experiment was carried out at The Svedberg Laboratory (TSLi in Uppsala in February 2005 using a setup 
called MEDLEY. The MEDLEY setup is designed to measure light ion reactions like (n,p), (n,d), (n,t), (n,'^He) 
and (n,a). Besides calcium, MEDLEY has also performed cross section measurements at 96 MeV on other 
nuclei of biological or technical relevance like silicon [22], oxygen ||231 . carbon |24J, iron, lead and uranium 
ES j . This chapter contains a brief description of TSL the neutron beam, and a more in depth description of the 
MEDLEY spectrometer setup. The last two sections are devoted to the electronics and data acquisition system 
involved. 



3.1 The Svedberg Laboratory 

The protons for the primary beam was produced by the Gustaf Werner isochronous cyclotron that can accelerate 
protons up to energies of 180 MeV and heavy ions with charge Q and mass number A up to 192 Q'^/A MeV. It 
consists of two D-shaped electrodes with a perpendicular magnetic field that makes the particles go in a circular 
pattern, and a high-voltage field between the electrodes that make the particles accelerate in the gap between 



the electrodes. By using a radio-frequency (RFi alternating electric field the particle will get a spiral trajectory 
outwards until extracted. The Gustaf Werner cyclotron is operated with a |RF| of 58 ns for 98 MeV protons. 
However, for protons of energies higher than 100 MeV, as well as '^He at their highest energies, relativistic 
effects need to be compensated for and the cyclotron switches to a frequency modulated synchrocyclotron 



mode where the driving RF is varied. This switching between modes is something that makes the Gustaf 
Werner cyclotron unique. 

From the cyclotron the proton beam is transported to the Blue Hall, as seen in the drawing in figure |3.1| 
Here it hits a highly enriched target of ^Li, where it induces the nuclear reaction ^Li(p,n)^Be with a Q value of 
1.64 MeV 1 3 1, resulting in a quasi-monoenergetic neutron beam in the range of 20-180 MeV. There are different 
thicknesses of the lithium target available, ranging from 2 to 24 mm ll26l . in this experiment the thickness 



used was 8 mm. The remaining proton beam is bent away using a bending magnet, actually recycled from the 



former LISA spectrometer, into a well-shielded beam dump (BD i and integrated in a Faraday cup. The charge 
deposited in the |BD| is used as one of three independent monitors for the neutron flux. 

The actual beam is defined by a collimator, consisting of 100 cm long iron cylinders of different diameter, 
embedded in a wall of concrete. These cylinders can give a beam of, more or less, any shape and diameter up 
to 30 cm, even though the standard forms available are circular beams in steps of about 5 cm. Right after the 
collimator is the experimental area with the MEDLEY setup, that is used for experiments concerning neutron- 
induced charged particle production and will be described in more detail in section |3.2| It is centered at a 
distance of 3.74 m from the lithium production target. After MEDLEY there is another setup called SCANDAL 
IIZT I used for experiments on elastic neutron scattering. Running several experiments at once like this is possible 
since almost the entire neutron beam passes the first setup without interacting. Two more independent neutron 



monitors are available namely the ionization chamber monitor ( ICM i that uses the ionization of the gas inside the 



detector and via an electric field collects the ion pairs created |28 1, and the thin-film breakdown counter ( TFBC| l 



. The exit area from the collimator as well as the |lCM| and the |TFBC| can be seen in figure [372| However, 



in this experiment there was a problem with the ICM so that it could not be used as an absolute monitor of the 



neutron flux, but instead as a relative monitor calibrated to the |TFBC| 

3.2 MEDLEY spectrometer 



The MEDLEY spectrometer, as seen in figure 3.3 consists of a 24 cm high cylindrical vacuum chamber with 
90 cm diameter. Inside this chamber there are eight telescopes mounted at even intervals covering angles from 
20 to 160 degrees, four on each side, and can be changed by rotating the table inside the chamber The distance 
from the target to the different telescopes can also be varied. 

To get a clean positioning of the target within the beam without getting too much background signal from 
the target holder, a small aluminum frame is used that is well off the beam. The target itself is strung up in 
the frame using thin wires. This setup can be seen in figure |3.4| There are three of these frames so one can 
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easily switch between targets without breaking the vacuum. For calibration purposes an ^^^Am source is also 
available, though this was not used in this experiment. 

3.3 Detector telescopes 



Each of the telescopes, as can be seen in figure 3.5 in the MEDLEY spectrometer consists of three individual 



detectors. First there are two silicon detectors, one thin of about 50 - 60 pm and one thick of about 400 - 500 pm 



thickness respectively. Finally there is a 5 cm long crystal of thallium activated cesium iodide (Csl(Tl) i. These 
detectors are often referred to as the AEi, the AE2, and the E detector, or sometimes A, B, and C detector. 
The detectors are mounted in an aluminum housing. Since the configuration of the detectors in this experiment 
was not the usual for MEDLEY it is listed in table lA.ll 

There also exists a possibility to use plastic scintillators mounted as active anti-coincidence collimators to 
discard the signals from particles that did not pass straight into the detector. Attempts have actually been made 
in earlier experiments, but due to problems with them they were not used in this experiment. 

3.3.1 Silicon detectors 

The two most common semiconductor detectors are silicon and germanium detectors, and are in principle a 
large diode in reverse bias. The primary way of energy loss is via Coulomb interaction in the material where the 
atoms are either ionized or excited. This will create electron-hole pairs that can be collected by two electrodes 
and measured as a small current pulse when a particle is passing through. 

Another feature that is fundamental for this experiment is that the amount of energy lost per unit length is 
not constant. Actually it decreases with increasing energy. In figure [33] the stopping power and the distance a 
proton can travel before it is completely stopped, the so called range, is plotted for the range of proton energies 
relevant in this experiment. This property will be used for particle identification in the analysis. 

These detectors are standard fully depleted surface-barrier detectors from ORTEC 1311 were the silicon 
wafer is glued to a metal ring. One side is covered with a thin layer of gold and the other with aluminum (3T\. 
The AEi detector is used as a transmission detector, that is the particles pass through the detector leaving only 
a fraction of its energy. This is also the case for the AE2 detector at energies above 10 MeV, but for low energy 
events the proton is stopped in this detector. Protons that stop in AEi are below the cutoff and are not included 



in the analysis. For an illustration of the different kinds of events, see figure 3.5 



3.3.2 CsI(Tl) detectors 

An inorganic scintillator is used as E detector in this setup. Here the incoming radiation creates electron-hole 
pairs in the crystal by moving an electron from the valence band to the conduction band on its way through the 
crystal. In the crystal, some impurities are often added. These impurities will have lower ionization energy than 
the crystal, which is why the holes will drift here and ionize the impurities. The ionized impurities will then 
pick up an electron from the conduction band and decay to the ground state and emit a scintillation photon f28l. 
A quite popular inorganic scintillator is the cesium iodide, available with both thallium and sodium acti- 



vation. The high density of the CsI(Tl) gives it large stopping power and it also has a really high light yield 
of about 65000 photons/Me V |28|. The wavelength of the emission light, 540 nm, is quite bad suited for a 
photomultiplier tube so instead a Hamamatsu photodiode is used |31 1. To suite this, the last 2 cm of the 4 cm 



diameter cylinder is tapered off to 1.8 cm to fit the surface of the diode. Another useful feature of the CsI(Tl) 
worth mentioning is that it has a variable decay time for different kinds of particles [28 ] that makes it suitable 
to use with pulse shape analysis. 

3.4 Readout and electronics 

Both the AEi and the AE2 detector can generate a trigger signal, in order to allow triggers from both low- and 
high-energy particles. The signal from the detector is split up in two branches, a timing branch and an energy 



branch. A complete scheme is found in figure 3.7 
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A signal from a detector is first amplified and then, if the pulse height is enough, converted to a logical 



signal using a constant fraction discriminator (CFDi. The timing signal is split up into two branches. One of 



them is connected via a fan-in/fan-out ( FIFO i logic to the other detector and telescopes as a master event. The 



other branch is delayed and sent as a stop signal for the timing channel in the time-to-digital converter ( TDC[l. If 
the computer is busy processing earlier event the current event is discarded, otherwise it, through a logic FIFO| 



generates an energy gate signal for the analog-to-digital converter ( ADC I and charge-to-digital converter ( QDC i 



units, a start signal for the TDC and finally a trigger for the data acquisition. 

The energy signal is recorded from all three detectors and all eight telescopes as soon as one detector 



triggers. It is sent into a ADC with the previously named master signal as gate, after amplifying and shaping. 
From the cyclotron one also gets a |RF| signal that indicates a pulse from the neutron beam. The|RF|is split up 



the time-of-flight ( TOF i channels in the 



TDC 



through a|FIFO|in to thre e branches. Tw o of the |RF| signals, where one is delayed, are used as a stop signal for 

signal as a veto for the master 



It is also possible to use this 



RF 



signal if the cyclotron is switched off. 



3.5 Data acquisition 

For the data acquisition a system called SVEDAQ is used, and it is summarized here from ll34l . The system is 
taken from the EUROGAM experiment and has been modified to be suitable for use at the |TSL| It is divided 
into three main blocks, as can be seen in figure 3.7 connected via two independent networks. The three blocks 



are the event builder, the tape server and the control and monitoring workstation. 

The event builder is the block that reads out the data from the lADCl units in the VME crate as well as the 
TDC units and scalers in the CAMAC crate. The event builder itself is also located in the VME. It mainly 



consists of a Motorola CPU that runs the readout program and sends packages to the data network when the 
readout buffers are full. It also sends the computer busy signals for vetoing of the master event and the 100 Hz 
clock so that events occurring during processing of other events are discarded. 

The tape server is an optional part of the acquisition system. Of course if one wants to use the data after the 
experiment one needs to write it to the hard drive, but if one wants to conduct different kinds of test runs it is 
not needed by the event builder. 

To control all this a SUN workstation is used. It is via this workstation that the system is started and 
stopped, but since the three blocks are more or less independent of each other the event builder and tape server 
will continue to write data from the CAMAC and VME even if the workstation is completely shut down. This 



workstation is connected to the local area network at TSL and can be remotely controlled from there. But it is 



also connected to the data network via a sort-spy daemon process so it can listen to the traffic without formally 
being a part of the network with its own IP number This makes it possible to display a large number of one- 
and two-dimensional spectra from the CAMAC signals and conduct some on-line sorting and analysis. 



'it may sound backwards to use the registration of a particle as a start signal and the next beam pulse as stop for the |TDc| but this is 
simply a matter of not overloading the electronics with unnecessary start signals. The reason to store the |RF| two times where one is delayed 
may not be intuitively clear either, but this is explained in more detail in section|4.3.2| and also in 1331 . 
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Figure 3.2: In the left panel, the exit of the collimator with the larger cylindrical \ICM\ and the smaller rectangular \TFBC\ in 
front. The right panel shows the neutron spectrum emerging from the collimator h301 . 




Figure 3.3: The MEDLEY spectrometer. The left panel shows a schematic sketch of the scattering chamber and in the right 
panel is a photography of a similar setup. 




Figure 3.4: One of the aluminum frames and the CH2 target stinged up inside the MEDLEY spectrometer. 
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Figure 3.6: Stopping power (solid line) and range (dashed line) of protons in silicon calculated using SRIM. 
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Figure 3.7: Scheme of the MEDLEY electronics as of week 9, 2004 



4 Data analysis 



The analysis of the data is split up into several parts. Since the SVEDAQ system writes the events as binary 
data, the raw files first have to be properly decoded. After the decoding, the information needs to be calibrated 
so that it corresponds to correct energy and timing information, before the particle selection procedure can take 
place, where events corresponding to the right particle type and incoming neutron energy are selected. Finally 
the relevant corrections for background and other effects can be carried out before extracting the final cross 
section. 

During the analysis, three different data sets were analyzed simultaneously. Two of them were the cal- 
cium data and the background data, while the third was a reference data set of polyethene for calibration and 
normalization purposes. The complete set of runs is listed in tables [A!2] and |A.3| 



4.1 ROOT 

The main tool used for the analysis is ROOT. ROOT is an object oriented data analysis framework initiated 
at CERN by Rene Brun and Eons Rademakerj*"] The framework is freely distributed as open source, where 
everyone is free to further distribute and modify it. This is something that has contributed, and still is contribut- 
ing, to its rapid growth. It contains a large amount of packages concerning different part of data analysis like 
histogramming, curve fitting, minimization, statistics tools and much more fSSl. 

The framework of ROOT is closely bound to a C++ interpreter called CINT that is written by Masahuru 
Goto. This means that one actually can use C-n- as a scripting language for rapid prototype development of 
programs. Then the same code can be compiled to take advantage of the fast running of a machine language 
executable, unlike if one had used a normal scripting language like Ruby, Perl or Python for the prototyping. 



4.2 Decoding 

The raw data from SVEDAQ are written into binary data files in the EUROGAM standards with a 24 byte 
header followed by event-by-event records |34|. A typical example of a data sequence can be found in figure 
The events are recorded as a series of short^^ with the first one being the number of bytes that are used 
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to register one event. Between the first short and the telescope events there are 20 scaler channels with numbers 
1-13 containing information about the telescope triggers, the 100 Hz clocks and the neutron monitors. The last 
seven channels are left empty. Next is the telescope events, each of them written in five shorts corresponding 
to three energy signals, one from each detector, and timing signals from the two silicon detectors. Finally the 
two RF signals are written before the event is ended. To translate the SVEDAQ files into ROOT files a small 



computer code called SV2ROOT 11361 was used 



4.3 Calibration 



4.3.1 Energy calibration 



Since the data from the runfile are given in the unit of channel numbers from the ADC an energy calibration is 
needed to convert these to incoming energies. The silicon detectors are assumed to have a linear response to the 
energy, that is the correlation looks like 

E^kx + m (4.1) 
where E is the energy output, x is the channel number, and k and m the constants to be fitted. 



Unfortunately the Csl(Tl) is not that simple since it has a non-linear response function between light output, 
L, and incoming energy, E, which makes the calibration more complicated. The light output is assumed to be 
given by a three-parameter formula that follows 

/ X 2i \E + a2Az'^' 

L — ao + ai [ E — a2Az In 



(4.2) 



'"These two have also been a part of creating other popular data analysis and simulation tools like PAW, PIAF, and GEANT. 
' ' One short is 2 bytes 
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Tl 






T2 


T7 




T8 






A£;i 


AE2 


AEi 


AE2 


AEi 


AE2 


AEi 


AE2 


Thickness [|jm] 


64.9 


549 


60.5 


538 


52.9 


424 


61.6 


550 


p [MeV] 


2.42 


8.62 


2.32 


8.51 


2.34 


8.63 


2.13 


7.41 


d [MeV] 


3.12 


11.5 


2.99 


11.4 


3.02 


11.5 


2.74 


9.87 


t [MeV] 


3.60 


13.5 


3.43 


13.4 


3.47 


13.6 


3.14 


11.6 


3He [MeV] 


8.60 


30.6 


8.23 


30.2 


8.32 


30.6 


7.56 


26.3 


^He [MeV] 


9.56 


34.4 


9.14 


34.0 


9.24 


34.5 


8.38 


29.6 


^Li [MeV] 


18.1 


65.3 


17.3 


64.6 


17.5 


65.4 


15.9 


56.3 


''Li [MeV] 


19.1 


69.7 


18.3 


68.9 


18.5 


69.8 


16.7 


60.0 



Table 4.1: Thickness and punch through energies for the different silicon detectors used. 



where is the fitting constants, and A and z is the mass number and charge of the incident particle |'3T1. 
However, to conduct the calibration we need the inverse of expression (4.2 1. Since this inverse is analytically 
complicated to calculate, the approximation according to Dangtip et al. [31] 



E 



+ hL + c{hLf 



(4.3) 



is used for the hydrogen isotopes. The parameter c depends only on particle type and is found by Tippawan lfT2l 
to be c = 0.0032 for protons, which leaves only a and h as constants to be fitted. 

Here the calibration of the silicon detectors is done by calculating the incoming energies needed to break 
through each of the silicon detectors - the punch through energies. These are calculated using a program called 
The Stopping and Range of Ions in Matte r (jSRIM 1, based on the work of Ziegler et al. OTl . The punch through 
values for Telescope 1 ( |T1| |, Telescope 2 ( |T2| i, Telescope 7 ( |T7[ ) and Telescope 8 ( |T8| l are found in table |4.1| 

If one plots the AEi versus AE2, these punch though points manifest themselves as turning points in the 
particle bands as can be seen in the left panel of figure 4.1 It is now quite straightforward to compare the 



channel numbers of these areas to the values in table 4.1 and get the coefficients to (4.1 1 




» e 10 12 14 16 
AE, [M«V1 



RUN015 




40 60 
E(MeVl 



Figure 4.1: Calibrated AE1-AE2 and AE2-E plots for \Tl\ The red lines correspond to tabulated energy loss values for 
protons, deutrons, tritons, "^He and alpha particles in ascending order The apparent mismatch for some lines in the right 
panel is due to the non-linear effects in the \Csl(Tl)\ 
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But once again, the CsI(Tl) poses a bit of a problem. As can be seen in the right panel of figure 4. 1 there 



are no turning points available for calibration. Instead, at least in forward angles, one can use the quite strong 
peak of elastic p(n,p) scattering in the CH2 runs as one calibration point. In the lower part of the E scale the 
proton band is quite bent, so here a couple of points are selected and the energy deposited in the E detector is 
calculated from the energy lost in AE2- Numerical values obtained for the calibration constants are hsted in 
table lA4l 



4.3.2 Time calibration 

For the time calibration a couple of special runs devoted were made. During these runs cables of different 



lengths, giving different delay times, were put between the RF signal and the TDC displacing the peak in the 



recorded RF signal. The result from these runs can be fitted with a linear calibration formula 



When each of the |RF| signals are time calibrated they need to be combined into a single time line. The two 
RF| signals can be plotted against each other in figure |4.2[ where four areas can be identified, schematically 
illustrated in figure |4.3| Events in area 1 and area 3 are events that occur during one of the |RF| signal pulses, 
and hence get start and stop signals simultaneously. So in these areas we only get signal from one of the|RF| In 



area 2 and area 4 there is a signal from both RF 1 and RF 2 



To translate the two RF signals into a single time scale, one of the |RF| signals is chosen as base, in this case 
|RF| 1 . In the areas where no information exists for the chosen base |RF| signal, the information from the other |RF| 
signal is transformed to the chosen base according to the algorithm in 1331 . 



RUN030 




RF2 [channel #] 



Figure 4.2: RF 1 versus RF 2. Each point in the 8- shaped trajectory correspond to a time between detection and production 
of the next burst. One lap around is 58 ns. 



4.4 Particle identification 



As mentioned earlier, the bands that are shown in figure [4~T| correspond to different kinds of charged particles. 
This makes the particle identification procedure quite straightforward by taking some rough cuts around the 
structures of the desired particle type. These cuts should be very generous, since it is always easier to remove 
unwanted events than to recreate wanted one ^and is shown for ' 
event must be in both of these to be selected. 



Tl 



in figure 



4.4 



The criterion here is that an 



'^At least concerning background events. Misidentified events from, for example, deuterons might still be tricky to get rid of if the cuts 
are too generous. 
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Figure 4.3: Scematic illustration of the two \RF\ signals and the different possible events that can occur h33'i . The |i^F| is 
measured as the from the trigger to the next signal from the cyclotron. 
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Figure 4.4: First graphical selection cuts for protons in \Tl\ In the left panel the cut for the AE1-AE2 plot is shown and in 
the right panel the AE2-E plot is shown. 



When these first rough cuts are done, the particle bands are straightened for a more refined selection of 
events. To do this, the energy recorded in the E detector is used together with |SR1M] to calculate the corre- 
sponding energies in the AEi and AE2 detectors. By subtracting this value from the energy actually recorded 
the particle bands are straightened out around zero. The A£^i(exp) — A_Bi(tab) selection is mainly used to 
remove noise in the low energy region in the AE1-AE2 plot, and the Ai?2(exp) — Ai?2(tab) selection is used 
to get a cleaner separation between protons and deuterons. These two are shown in figure [43] 



4.5 Time-of-Flight measurements 

When only protons are selected, the next restriction is to only select protons that are produced by neutrons of 
the correct energy. Since the neutron beam is not completely monoenergetic, but contains about 60 % lower 



energy neutrons [31 1, the contribution from these has to be removed. To do this the TOF information obtained 
in section f 



calculated as 



4.3.2 is used. The TOF t, for a particle with energy E^ and velocity v to travel a distance s can be 

t=-. (4.4) 



The velocity of the particle is given in its rest mass, ttiq, and its momentum, p, as 

P 



(4.5) 
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Figure 4.5: Refined selection cuts for protons in \Tl\ In the upper left panel the cut for the AEi{exp) — AEi{tab) plot is 
shown and in the upper right panel the AE2{exp) — AE2{tab) plot is shown. The lower panels contains the final proton 
cuts as in figure \4A\ for \Tl\ together with the tabulated curve from \SRIM\ 



where the p is given by 



1 



Plugging (4.6 1 into (4.5 1 into (4.4i and simpUfying yields 

s / {m,oc'^ + Ek)^ 



t 



(4.6) 



(4.7a) 



+ ^ 



= ^ jm^ + , 2 "^0 



(4.7b) 



where Eqs. (4.7ai and (4.7b i are equivalent 
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But the time scale mentioned in section 4.3.2 includes both the neutrons TOF as well as the charged particles 



TOF from the target to the telescope. As can be calculated from (4.7ai, if the particle is a 94.5 MeV neutron 



traveling a distance of 3.74 m the TOF would be about 29.9 ns. A 5 MeV proton traveling 26.83 cm has a TOF 



of 8.71 ns while a 90 MeV proton would have a TOF of only 2.19 ns. This makes quite a large difference 



especially at low energies, and is compensated for by subtracting the charged particle |TOF| calculated from the 



energy of AE'i, Ai?2 and E added together The final TOF cuts are illustrated in figure [4^6 
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Figure 4.6: Neutron \TOF\ versus proton energy for calcium in the left panel, and CH2 in the right panel, together with the 
\TOF\ cut. The the low energy tail in the neutron spectra is clearly visible in the CH2 plot. The extra \TOF\ cut at high energies 
is motivated in section\6^ 



To choose a suitable sized cut, the elastic p(n,p) scattering is chosen to define the width. This peak has a 
full width at half maximum (jFWHMj) of 6.6 ns, quite in agreement with Dangtip et al. [31 1 where the |FWHM| 
is measured to be about 6-7 ns, and one of the large contributions is the finite time of a beam pulse from the 
cyclotron. The with of the cut is chosen as two standard deviations, and is a compromise between statistics, 
minimizing mismatching errors, and a high lower limit of the accepted neutron spectrum. For more details, see 
section l64l 



4.6 Corrections 



For each run, the live time of the data acquisition system is corrected for. As mentioned in section [J?] events 
that occur while an earlier event is treated are discarded. The amount of discarded events is not constant, as seen 



in figure 4.7 but is assumed to be. The ratio of discarded events is obtained by comparing the 100 Hz scaler and 



the 100 Hz computer busy scaler. The finite efficiency of the CsI(Tl) is also corrected for This has a non-linear 
dependence of energy deposited, and this is simulated by both MCNPX and GEANT |38|. The results of these 
simulations are found in figure [4~8| Since they are in good agreement, one of them is choosen, in this case the 
GEANT results. 

During the entire event-by-event analysis, the data sets from the runs containing an empty frame were treated 
in parallel to the calcium and the polyethene data. When switching to analysis of histograms instead of events, 
the data from the background runs kan be subtracted from the calcium and polyethene data. In that way one can 
get pure signal histograms. A comparision between the calcium, polyethene and background data can be found 
in figure [49| The background in MEDLEY has been analyzed thoroughly by Tippawan 1 12 1 for the old facility. 
Its main components are particles produced by produced neutrons reacting in the beam pipe and the reaction 



chamber, and the CsI(Tl) and the telescope housing acting as an active target for the beam halo. 
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Figure 4.7: Typical live time of a run. The instabilities Figure 4.8: Efficiency of the \CsI(Tl)\ detector The 
in the beginning are due to statistics. dashed line is a simulation with MCNPX and the solid 

line is a simulation with GEANT [38 ]. 



Ion Energy [MeV] 


dE/dx [keV/^m] 


Range [^m] 


Straggling [^m] 


3.25 


11.854 


164.35 


7.51 


3.50 


11.247 


185.84 


8.35 


3.75 


10.708 


208.46 


9.19 


4.00 


10.225 


232.18 


10.06 


4.50 


9.3934 


282.85 


12.85 


5.00 


8.7004 


337.80 


15.54 


5.50 


8.1137 


396.92 


18.20 


6.00 


7.6086 


460.15 


20.86 


6.50 


7.1698 


527.41 


23.55 



Table 4.2: Stopping of protons in calcium 



4.6.1 Target thickness 

To achieve a count rate that is acceptable, some demands are put on the target. Since the neutron beam is 
quite low in intensity and only a small fraction of the neutrons interacts the target has to have some finite 
thickness. For high-energy protons this is of negligible importance. For low energy protons, on the other hand. 



the energy loss can be quite significant. In table 4.2 SRIM calculations for low energies are shown to illustrate 



the problem. The target used in these measurements was about 230 /im thick. These effects are compensated 
for through a method implemented in a computer code, TCORR, by Pomp and Tippawan ll39l and is based on 



iterative calculations of response functions. The result of the target correction procedure is found in figure 4.10 



4.7 Accepted neutron spectrum 

Even with the |TOF| selection in section [43] not all low energy neutrons were rejected. As can be seen in figure 



4.6 the finite width of the TOF cut results in events induced by neutrons down to around 60 MeV to be accepted. 
In the same figure one can also see the how the low energy p(n,p) tail wraps around the time scale and allows 
events by about 10-14 MeV neutrons inside the cut. To be able to subtract data originating from these events 
the accepted neutron spectrum is deduced from the CH2 data. To get a clean neutron spectrum, the contribution 
from carbon has to be subtracted from the CH2 data. For 95 MeV neutrons, the subtracted carbon cross section 
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Figure 4.9: In the top panels, the signal and background for calcium, left, and polyethene, right, at an angle of 20 degrees. 
The bottom panel shows the background for the four angles analyzed: 140 degrees, 160 degrees, 40 degrees and 20 degrees 
in increasing order. 



comes from existing experimental data f24l. This is not possible for lower energy neutrons, so their contribution 
to the carbon spectrum is taken from tabulated data in [40 1. 

To obtain a p(n,p ) spectrum, the peak at 83 MeV is normalized to a cross secti on ob tained from the partial 



wave analysis (PWAi solution SP05 of |4T1 via a procedure described in appendix e[ 



In this first normaliza- 
tion, the fact that also carbon contributes slightly to the peak is ignored, but the peak is re-normalized to the 
cross section from |41 1 after the subtraction of experimental carbon data. 

Now the remaining CH2 data are divided into bins corresponding to 10 MeV neutron energy each, translated 
to elastically scattered proton energies as seen in table A. 5 The bin right below the peak is integrated and the 



number of events in the bin is corrected for the variation in p(n,p) cross section, also listed in table A. 5 The 



'^One should note that this is not the same cross section as the one used for later normalization, but this should pose no problem since 
this is only a relative measurement. 
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Figure 4.10: Target correction results for calcium in the left panel and polyethene in the right panel. Red circles represent 
the original data, while blue triangles show the corrected data, re-binned to I MeV for comparison. After three iterations 
the Kolmogorov test gave 100% probability for the current iteration compared to the previous one. 



cross section is assumed constant within the bin and equal to the cross section in the center of the bin. When 
doing this an assumption is made that the high Q value of the C(n,p) reactiorp]makes the contribution from the 
carbon to the bin content negligible. Now, when one have the fraction of accepted neutrons in that bin compared 
to the full energy peak, one can subtract the corresponding carbon cross section obtained from |40|. In parallel, 
the obtained bin content relative to the peak is used to build up an accepted neutron spectrum. This process is 
then iterated downwards in energy for each of the bins until a pure p(n,p) spectrum and a corresponding neutron 



spectrum is available, and the resulting spectrum is found in figure 4.1 1 and in table 4.3 
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Figure 4.11: The deduced pure np spectrum in the left panel. Grey triangles are polyethene data and black circles are 
carbon data from the procedure in section^~7\ In the right panel is, the fraction of accepted neutrons per bin is shown. 



"The Q value for the C(n,p) reaction is 12.588 MeV gj. 



26 



Bin [MeV] 


Contribution [%] 


5-15 


3.2 


15-25 





25-35 





35-45 





45-55 





55-65 


6.9 


65-75 


10 


75-85 


17 


85-100 


61 



Table 4.3: Accepted neutron spectrum, as seen infigure\4.11\ 



4.8 Normalization 

To obtain a value for the cross section all data are normalized to elastic p(n,p) scattering at 20 degrees. This 
since the p(n,p) scattering cross section is well known and the peak is clear at this angle. The peak together with 



a gaussian fit is seen in figure 4.12 




Figure 4.12: The elastic np peak fitted with a gaussian. The data are CH2 with carbon H24\l subtracted. 



The FWHM of the peak is found to be 6.4 MeV, in agreement with the FWHM in Dangtip et al 
the solid angle of the telescope as the main contribution. The cross section, a, for each bin, i 



is now calculated from 



31), with 
in telescope TX 

(4.8) 



Ncs.A,TX Nil AfcH2 'TlCa ^TX ^Cr 

where is the number of counts, m is the target mass, AI is t he mole cular mass, $ is the relative neutron 



Tl 



and 



T8 are the same, as seen in table 



A.l 



the 



flu)|^and fl is the solid angle. Since the target distances for 
solid angles for these telescopes are assumed to be equal. The other six telescopes, as also seen in table |A.l 
however needs to be corrected with a factor of 0.546, from a simple geometrical point of view. For a more 
in-deep analysis of the solid angle in MEDLEY, see [42|. 



''The nonnalization to neutron flux has actually already been done, simultaneous to the life time correction, using the |lCM| monitor. But 
the factor should still be included in Eq. j4.8t for completeness. 



27 



Since the polyethylene is the reference target that is regularly used in MEDLEY experiments the weight, 
and number of protons, has been measured with high precision to be 461.55 ± 0.01 mg and 3.963 • 10^^ protons 
ll42l . The weight of the calcium target used was measured to be 235.9 ±0.1 mg before the runs and 237.3 ± 0.1 
after the runs, which averages to 236.6 mg. 



The cross section, ct™, is taken from Rahm et al. B3I where the p(n,p) cross section at a center of mass ( CM i 
angle of 139.0 degrees and at 96 MeV is measured to be 7.735 mb/sr. In order to extrapolate this to 94.5 MeV 
CM angle of 139.1 degrees, data from PWA solution SP05 are used. From (41] the cross section for 94.5 



at a 



MeV and 139.1 degrees is 8.041 mb/sr, while the cross section for 96 MeV and 139.0 degrees is 7.928 mb/sr. 
This gives an increase in cross section of about 1.43 %. Via the procedure in section [B] this is translated into 
crJ?'' = 30.53 mb/sr 



28 



5 Results 



The results given here are those of the accepted neutron spectrum obtained in section 4.7 and listed in table 4.3 
resulting in a mean value for the neutron energy of 79.7 MeV. The total cross section for this spectrum is found 
in table |A.12| with a cutoff energy of 2.5 MeV. 



5.1 Double-differential cross sections 

Double-differential cross sections at four different angles have been analyzed. The result is found in figure [5T| 
and tables [Ar6l|A.9| which also include corrected data introduced in section [F.l. 3 1 



9=20 degrees 




10 20 30 40 50 60 70 80 90 
Proton energy [MeV] 



9=40 degrees 




10 20 30 40 50 60 70 80 90 
Proton energy [MeV] 



9=140 degrees 



9=160 degrees 




10 20 30 40 50 60 70 80 90 
Proton energy [MeV] 




10 20 30 40 50 60 70 80 90 
Proton energy [MeV] 



Figure 5.1: Experimental double-differential cross sections, filled circles, for the neutron spectrum in table pTj] in the 
Ca(n,px) reaction. The solid histogram is TALYS calculations for this neutron spectrum. Cutoff energy is 2.5 MeV. 
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5.2 Single-differential cross sections 



5.2.1 Angular-differential cross section 

To get angular-differential cross sections, the double-differential spectra are summed over and weighted for the 



different bin widths. The result is found in table A. 11 In figure 5.2 the angular-differential cross section for 
bins corresponding to compound, pre-equilibrium and direct reactions is found. 



Compound reactions, 3-11 MeV [ 
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Figure 5.2: Experimental angular-differential cross sections for three different bins and the neutron spectra in table \4.3\ in 
the Ca(n,px) reaction. The chosen bins are assumed to represent the three specific reaction types. The solid line results from 
TALYS and the dotted line is a fit to the data according to Eq. \2.11\. Cutoff energy is 2.5 MeV. 



5.2.2 Energy-differential cross section 



To get energy-differential cross sections, each of the energy bins in the double-differential spectra are fitted with 
Eq. ( |2.1 Using the results from this fit, the data are interpolated to the angles 60, 80, 100 and 120 degrees, 
and also extrapolated to 5 and 175 degrees. These data points are summed over and weighted for the different 



bin as well as their angular contribution. The result is found in figure 5.3 and in table A. 10 
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Figure 5.3: Experimental energy-differential cross sections, filled circles, for the neutron spectrum in table \4.3\ in the 
Ca(n,px) reaction. The solid histogram is TALYS calculations for this neutron spectrum. Cutoff energy is 2.5 MeV. 
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6 Discussion 



6.1 Comparison with theory 

To be able to compare the result with theoretical calculations the contribution from the low-energy neutrons 
need to be subtracted. To do this three different methods, listed in sections 6.1.1|6.0 are presented. The 
difference in total cross section between these methods is shown in table 16. II 

6.1.1 Method 1: Rescaling of experimental data 

In the first method for subtraction of low-energy neutron data, the energy spectrum for proton emission is 
assumed to be independent of the incoming neutron energy. The experimental cross section is weighted with 



the amount of low energy neutrons in the accepted spectrum from table 4.3 and rescaled to the peak in the 



neutron spectrum. The result is found in figure 6.6 



The main advantage of this method is that one gets rid of the model dependencies in the TALYS calculations, 
and the data used are pure experimental. This means that the result at low energies, up to about 40 MeV, should 
be quite accurate. The main disadvantage on the other hand is the obvious fact that the cross section spectrum 
for 94 MeV neutrons and lower energy neutrons differ a lot at high energies. This is especially relevant at 
energies above the low energy bin under consideration. 

6.1.2 Method 2: Subtraction of TALYS data 

Another way to correct for the low-energy neutron contribution is by using data from TALYS. For the accepted 
neutron spectrum in table |4.3| TALYS is used to calculate cross sections for each bin. The bins between 15 MeV 
and 55 MeV are left out, since the contribution in these bins probably derive from statistics. The result is found 



in figure 6.7 



The advantage of using this method is that it makes the comparison between experiment and theory straight- 



forward. The difference one can see in figure 6.7 between experiment and theory is directly related to the TALYS 
calculations. So for code verification purposes, this method works fine. But for obtaining experimental data, 
this method is not so good, since the resulting experimental data points will have a heavy model dependency on 
calculations that are known to be incorrect. 

6.1.3 Method 3: Subtraction of modified TALYS data 

A third possible method is to try and modify the TALYS data to get it more correct. From previous experi- 
ment, TALYS is known to heavily overestimate the evaporation cross section in backward angles, and slightly 
overestimate it at forward angles. It also appears to underestimate pre-equilibrium cross sections in forward 
angles 1221 l23l . Here, an assumption is made that the over/under-estimations in the TALYS calculations are 
truly systematic and do not vary with projectile energy. From this assumption one can divide the experimental 
data with the TALYS data in figure [5T| and get a correction histogram. This correction histogram is assumed to 
contain the true over/under-estimations in TALYS, and by multiplying the low-energy TALYS calculations with 



this before subtracting them, the truq^cross-section spectra are obtained. The result is found in figure 6.8 and, 
together with uncorrected data, in tables |A.6[|A.9| 

By using this method one can reduce the model dependencies compared to unmodified TALYS data, by 
adjusting the theory to the experiment. This method is the one that probably will give the result with smallest 
systematic errors for the low-energy neutron correction. However, it is first when one sees the results from 
similar experiments at different energies that one can conclude if the assumption made is too bold or not. 

6.2 Background and shielding 



As can be seen in figure 4.9 the background problem is quite heavy. In the silicon data runs from the old facility. 



the signal-to-background ratio was about 8 lfT2l . By moving the spectrometer closer to the lithium target, and at 



'*Well, true under the condition that the assumption is correct. 
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Method 


3^(^ = 20) 






1^(0=160) 


CTprod(n,p2;) 


1. Exp. Data 


155 ±2.92 


99.4 ± 1.54 


31.6 ±0.486 


32.1 ±0.700 


746 ± 1.97 


2. TALYS 


188 ±4.76 


107 ±2.51 


17.6 ±0.792 


18.6 ± 1.14 


679 ± 3.74 


3. Mod. TALYS 


175 ±6.69 


107 ±4.22 


33.7 ± 1.73 


35.1 ±2.04 


788 ± 5.83 


TALYS 


131 


106 


62.4 


62.4 


1070 


GNASH ilPI 


186 


145 






1300 



Table 6. 1 : Comparison of angular-differential and total cross sections with subtracted data. Differential cross sections 
are given in mb/sr and total cross sections in mb. The cutoff energy is 2.5 MeV, except for the GNASH calculations of 
o'pmd{n,px), where no cut-off is was applied. 



the same time reducing the shielding, the background increased dramatically. In later data sets with energies of 
180 MeV, the background problem were so large that it was virtually impossible to get acceptable results from 
the data |36|. To come to terms with the background problems recent GEANT simulations at St. Petersburg 
University show that the background for a 180 MeV neutron beam can be reduced by a factor 20-5Cp^ by 
exchanging the concrete wall to an iron wall. The background would probably be reduced even more for a 94 
MeV beam. Available iron blocks from the old CELSIUS ring have actually been used when this reconstruction 
was undertaken in January 2007. The first data with the new shielding will probably be taken during February 
2007. 



6.3 Special TOF cut 



Run 37 was itself unusable for the analysis, but still it may contain some interesting information on the phe- 
nomena. As can be seen in the right panel of figure 6.1 there exists a really strong shift in the TOF at high 
energies that does not seem to be constant, but disappears around 60 MeV. This effect is also visible in Tl as 



can be seen in the left panel of figure 6.1 This clear energy dependence makes a correction for this effect quite 



difficult. Therefore, the TOF cut instead has been expanded to include these effects. This had to be done for 



both Tl and T8 since the effect still was present, but weak, after the threshold adjustment in |T8| The reasons for 
this effect is unfortunately still unknown, but since it was weakened by threshold adjustmenj"^|it is suspected to 
be correlated to the pulse height in the AE2 detector This suspicion is further strengthened by the fact that an 
analysis of emitted deutrons shows no effect like this. 

One can also note the energy shift in the lower part of the plot, that does not seem to be dependent on the 



threshold adjustment. If the cause of these problems is because of instabilities in the CFD a possibility to get 
rid of this pro blem is to use a timer based on the leading edge technique instead. This will require some extra 
work with the TOF selections and mean large corrections^ since the timing of a leading edge discriminator is 



dependent on the pulse height. 



6.4 Mismatch between Tl and T8 



As one can see in figure 6.2 there are some mismatches between Tl and T8 The low energy mismatch has been 



carefully examined, unfortunately without solving the problem. The intermediate energy mismatch seen in the 



calcium data in figure 



upper part of the TOF 



6.2| was however found to depend on the |TOF| cut. The mismatch was reduced when the 



cut in figure 4.6 was increased. But increasing the TOF cut also increases the amount of 



accepted low energy neutrons that need to be subtracted from the data. So the chosen TOF cut of two standard 
deviations from the peak, is a compromise between the error from the telescope matching and the error from 
the low-energy neutron correction. 



"Even up to a factor 100. 

'^Before njn 39 the AE2 threshold in T8 was changed from 450 mV to 315 mV, and the AE2 threshold in T7 was changed from 400 
mV to 325 mV. Before run 40 the AE2 threshold in T8 was further changed to 295 mV, and the AE2 threshold in T3 was changed from 
320 mV to 350 mV. 

"Of the order of 10 ns. 
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Figure 6.1: TOF shift in\Tl\in the left panel and run 37 in\T8\in the right panel. 
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Figure 6.2: Mismatch between [77] (filled circles) and \T8\ (filled triangles). The left panel shows the CH2 runs and the right 
panel shows the calcium runs. 



6.5 Reaction losses 

Due to nuclear interactions within the CsI(Tl) some protons may, for example, produce gamma radiation that 
will escape the CsI(Tl) undetected. Thus only a certain fraction of the incident protons is detected at its true 
energy, and the rest will form a reaction tail of apparently lower energy. This issue is currently under investi- 
gation, but so far the results show that for low energies the lost fraction is really small and correcting for these 



effects only gets important at higher energies, as shown in figure 4.8 
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Origin Uncertainty 

Target correction 1-10% 

Solid angle 1-5 % 

Beam monitoring 2-3 % 

Number of Ca 1% 



CsI(Tl) efficiency 1 % 

Particle identification 1 % 

Dead time < 0.1 % 
Absolute cross section 5 % 

Table 6.2: Systematic uncertainties 



6.6 Different approaches to angular distributions 



In this work, the angular distributions were obtained by fitting Eq. ( |2.1 1[ ) to the data. This is an approach 
coming from experiment, and trying to fit some angular shape to the data, that can be integrated to get the cross 
section. Another approach one can use is to go from theory and calculate the predicted shape, a, of the angular 
distributions and fitting Eq. (2.10i to the data under assumption that /msd = 1. The actual form of /msd, as 
calculated by TALYS, is found in figure [64| To calculate a, the parametrization used fSOl is 



a = 0.04 



1.8 • 10" 



Eieh' 



6.7-10 



_7 ( Eaet' 



(6.1) 



where Ei — min(ea', 130), ~ min(ea',41), ej = Ea + Sa and e^' — Ei, + St- Here Ea is the incoming 
neutron energy, Eb is the outgoing proton energy, and Sa,b is the separation energy between particle a, b and 
the compound nucleus. The separation energy, Sa,b, is calculated using the Myers and Swiatecki mass formula 
lll4l . As seen in figure [63] there exist some differences in angular distributions between the theoretical and the 
experimental approach, especially at higher energies, where the theory predicts the forward peaking be larger. 
The negative values on a at really high energies, that would imply backward peaking, are probably statistical 
artifacts from the background subtraction. 



As can be seen in figure 6.5 this parametrization gives a result closer to the model at intermediate energies 
but deviates more at high energies, while at low energies the difference between the methods are small. This is 
quite expected from the approximation of /msd = 1. An investigation of the parametrization with the calculated 
form of /msdj or even experimental determination of /msd. is something 1 think could be interesting in the future. 



6.7 Systematic uncertainties 

Throughout this work the error bars only represent the statistical uncertainties in the experiment. To give a 



complete picture, the systematic errors should be added to the result. These are listed in table 6.2 where some 
values are for this particular experiment, and some values are adapted from Tippawan et al. 1,23 J . 



6.8 Conclusions 



As presented in sections 6.1.1|6.13 and shown in figures 6.6 - 6.8 comparisons between three different methods 
have been made. No matter which method one chooses, the results all show the same trend as the previous 
experiments Il22ll23l . The evaporation peak is overestimated by TALYS and predicted to be almost isotropic, 
while in experiment it significantly decreases with angle. The GNASH calculations overestimate the peak even 
more. In the pre-equilibrium region, TALYS instead underestimates the cross section, but the predictions seem 
to be more accurate at larger angles. 



The problem in the TOE is an issue that needs to be investigated, especially since similar effects have been 



noted in other data sets as well |24|. But since one of the mismatches in section 6.4 was found to be dependent 



on the TOE cut, a solution to the TOE problem might shed some light on the mismatch. 
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Figure 6.5: Theoretical and experimental approaches to angular distributions and the differences in energy-differential 
cross sections for the two approaches, where experimental angular distributions are filled circles and calculated angu- 
lar distribution are filled triangles. The data used are the data with modified TALYS calculations subtracted. The solid 
line is prediction by TALYS, and the dashed line is prediction by GNASH h40'i . 



6.9 Outlook 

The data sets analyzed in this work was the first, and at this moment the only, data set available from the 



upgraded neutron facility |26] at TSL To further improve the calcium results, experimental data from the 
angles 60, 80, 100 and 120 degrees should be analyzed. The reactions Ca(n,d), Ca(n,t), Ca(n,'^He) and Ca(n,Q;) 
would also be interesting to have a look at. 

Apart from that, and as mentioned in section [3] lots of work has previously been carried out at 96 MeV 
II22II23II25I and there still exists a never-ending selection of isotopes to be measured. But since the nuclei under 
study in these works are especially interesting from a biological and technical point of view it is more interesting 
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to use the capabilities of the new neutron facility and ran the experiments at 180 MeV. By obtaining more data 
sets, one extends the relevant database for future applications. At the same time one makes it possible to do a 
more quantitative analysis on the nuclei and get a better systematic feedback for the development of TALYS. 
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Figure 6.6: Double-differential cross sections in the four top panels, and angle integrated energy-differential cross section 
in the bottom panel, with experimental data (method 1) subtracted for 94 MeV neutrons in the Ca(n,px) reaction. Solid lines 
are predictions by TALYS, and dashed lines are predictions by GNASH h40'S . Cutoff energy is 2.5 MeV. 
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Figure 6. 7: Double-dijferential cross sections in the four top panels, and angle integrated single-differential cross section 
in the bottom panel, with TALYS data (method 2) subtracted for 94 MeV neutrons in the Ca(n,px) reaction. Solid lines are 
predictions by TALYS, and dashed lines are predictions by GNASH 140J. Cutoff energy is 2.5 MeV. 
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Figure 6.8: Double-differential cross sections in the four top panels, and angle integrated energy-differential cross section 
in the bottom panel, with modified TALYS data (method 3) subtracted for 94 MeV neutrons in the Ca(n,px) reaction. Solid 
lines are predictions by TALYS, and dashed lines are predictions by GNASH h40'S . Cutoff energy is 2.5 MeV. 
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Appendices 

A Figures and tables 




Figure A.l: Contribution to total cross section from various reaction types, according to TALYS calculations. 



No 


Sil 


Si2 


Distance [cm] 


Angle [deg] 


Axi [/Ltm] 


/S.X2 [/um] 


1 


38-1 lOB 


35-150B 


26.83 


20 


64.9 


549 


2 


37-198D 


37-051B 


19.83 


40 


60.5 


538 


3 


38-1 lOD 


37-052A 


19.83 


60 


63.9 


533 


4 


37-198G 


37-052C 


19.83 


80 


61.7 


526 


5 


31-325C 


37-199A 


19.83 


100 


50.4 


439 


6 


31-325D 


37-199E 


19.83 


120 


50.1 


432 


7 


38-1 lOG 


39-lOOB 


19.83 


140 


61.6 


550 


8 


31-325H 


37-199G 


26.83 


160 


52.9 


424 



Table A.l: Detector setup 
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Run 


Target 


Beam 


Neutron fluence 


Lifetime 


Runtime 






[pA] 


[n/cm2 • IQS] 


[%] 


[h] 


004 


Junk 










005 


Empty 


1.2 


3.1 


84 


0:57 


006 


Empty 


1.1 


4.2 


74 


1:15 


007 


Empty 


1.2 


5.2 


72 


1:29 


008 


Empty 




5.5 


72 


1:35 


009 


Empty 








1:27 


010 


Empty 




5.9 


72 


1:36 


Oil 


Empty 


1.2 


3.5 


72 


0:57 


012 


Empty 




5.6 


72 


1:33 


013 


Empty 




6.9 


72 


1:55 


014 


CH2 


1.2 


6.5 


71 


1:45 


015 


CH2 




5.5 


71 


1:33 


016 


CH2 




4.3 


71 


1:09 


017 


Empty 


1.2 


4.7 


72 


1:22 


018 


Empty 


1.2 


4.7 


73 


1:26 


019 


Ca 


1.1 


1.6 


74 


0:27 


020 


Ca 




5.9 


74 


1:50 


021 


Ca 




7.1 


74 


2:02 


022 


Ca 


1.1 


7.1 


74 


2:11 


023 


Ca 




5.9 


74 


1:47 


024 


Ca 




5.6 


74 


1:37 


025 


Ca 




4.6 


73 


1:19 


026 


CH2 




5.1 


72 


1:28 


027 


CH2 




4.6 


72 


1:19 


028 


Ca 




7.1 


73 


2:01 


029 


Ca 




6.9 


73 


1:56 


030 


Ca 




5.5 


73 


2:00 


031 


Ca 




3.4 


73 


0:52 


032 


Ca 




5.7 


73 


1:37 


033 


Ca 




6.6 


73 


1:55 


034 


Ca 




7.0 


73 


2:01 


035 


Ca 




3.7 


74 


1:06 


036 


Ca 




1.6 


74 


0:35 



Table A.2: List of runs when the table was positioned at 180 degrees, that is with \T l\ in forward and \T8\ in backward direction. 
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Run 


Target 


Beam 


Neutron fluence 


Lifetime 


Runtime 






[pA] 


[n/cm2 • IQS] 


[%] 


[h] 


001 


Junk 










002 


Junk 


1.2 








003 


Junk 


1.2 








037 


CH2 


1.1 






1:05 


038 


Junk 










039 


CH2 


1.1 


0.96 


71 


0:15 


040 


CH2 




5.6 


71 


1:43 


041 


CH2 


1.1 


5.5 


71 


1:35 


042 


CH2 




1.8 


71 


0:30 


043 


CH2 




6.9 


71 


2:00 


044 


Empty 




5.9 


72 


1:44 


045 


Empty 




9.3 


72 


2:44 


046 


Empty 




7.2 


72 


2:07 


047 


Empty 




6.9 


72 


2:00 


048 


Empty 




7.0 


72 


2:01 


049 


Empty 




6.7 


72 


1:59 


050 


Empty 




6.8 


72 


2:00 


051 


Empty 




5.2 


72 


1:31 


052 


Ca 


1.1 


4.5 


71 


1:15 


053 


Ca 


1.1 


5.3 


72 


1:37 


054 


Ca 


1.1 


5.4 


72 


1:35 


055 


Ca 




5.7 


72 


1:45 


056 


Ca 




3.4 


72 


1:00 


057 


Ca 




4.7 


72 


1:25 


058 


Ca 




6.4 


71 


1:55 


059 


Ca 




6.9 


71 


2:02 


060 


Ca 




6.4 


71 


1:56 


061 


Ca 




6.6 


71 


2:01 


062 


Ca 




4.2 


72 


1:17 


063 


Ca 


1.1 


5.4 


72 


1:39 


064 


Ca 


1.1 


6.5 


72 


1:59 


065 


Ca 




6.5 


72 


1:59 


066 


Ca 




3.5 


72 


1:02 



Table A.3: List of runs when the table was positioned at degrees, that is with \T8\ in forward and \Tl\ in backward direction. 
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013C 


4264 


04FE 


0000 


0000 


4265 


05B2 


0000 


0000 


4266 


083E 


0000 


0000 


4267 


039B 


0000 


0000 


4268 


06B5 


0000 


0000 


4269 


06F4 


0000 


0000 


42 6A 


048B 


0000 


0000 


426B 


0325 


0000 


0000 


426C 


0A2C 


0000 


0000 


426D 


074E 


0000 


0000 


426E 


0C53 


0000 


0000 


426F 


OOOE 


0000 


0000 


4270 


02D9 


0000 


0000 


4271 


0000 


0000 


0000 


4272 


0000 


0000 


0000 


4273 


0000 


0000 


0000 


4274 


0000 


0000 


0000 


4275 


0000 


0000 


0000 


4276 


0000 


0000 


0000 


4277 


0000 


0000 


0000 


4501 


0066 


0065 


004F 


0000 


01F4 


4502 


004B 


02CA 


0137 


0000 


007E 


4503 


007A 


0068 


0040 


0000 


OOlC 


4504 


004B 


0034 


0043 


0000 


0206 


4505 


0064 


005A 


0036 


0000 


01F5 


4506 


0077 


0049 


003F 


0000 


OICF 


4507 


0062 


006D 


00D4 


0000 


0225 


4508 


0083 


0065 


004D 


0000 


01D3 


480A 


0060 


0068 


0067 


004D 


0069 


0047 


005A 


0062 


0000 


4C09 


0025 


0038 


004A 


004E 


0040 


0039 


002E 


002F 


OFDC 


OFDC 


OFDC 


OFDC 


0000 


4300 


0000 


0083 


018F 


FFFF 







Figure A.2: Typical sequence of hexadecimal raw data for an event. Group identifiers are marked in bold. Groups 42XX 
are the scalers, 45XX are the readout from the telescopes and the group 4300 0000 is the identifier for the \RF\ signals. The 
event is ended by FFFF. 





Tl 


T2 


T7 


T8 




0.00334 


0.00332 


0.00309 


0.00278 


mi 


-0.264 


-0.197 


-0.284 


-0.256 


k2 


0.00893 


0.00863 


0.00969 


0.00904 




-0.851 


-0.667 


-1.04 


-0.893 


a 


-2.09 


-0.734 


-1.57 


-1.49 


b 


0.0259 


0.0199 


0.0226 


0.0228 


c 


0.0032 


0.0032 


0.0032 


0.0032 



Table A.4: Calibration values, as used in Eq. \4.1\ and Eq. \4.3^ . 



Neutron bin [MeV] np bin [MeV] (^{n,^) [mb] 



5-15 


4-13 


292.60 


15-25 


13-22 


153.78 


25-35 


22-31 


100.51 


35-45 


31-40 


73.033 


45-55 


40-48 


57.336 


55-65 


48-57 


47.461 


65-75 


57-66 


40.841 


75-85 


66-78 


36.239 



Table A.5: Energy binning for subtraction of carbon data from the CH2 runs. 
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Bin 


Uncorrected 


Corrected 


Bin 


Uncorrected 


Corrected 


[MeV] 


[mb/sr/MeV] 


[MeV] 


[mb/sr/MeV] 


2.5-3 


0.0931 ±0.23 


0.0559 ±0.25 


68- 


-70 


2.67 ±0.30 


2.25 ±0.37 


3-4 


9.61 ±0.75 


5.66 ± 1.1 


70- 


-72 


2.13 ±0.29 


1.85 ±0.32 


4-5 


9.81 ±0.50 


5.90 ±0.86 


72- 


-74 


1.87 ±0.28 


1.67 ±0.30 


5-6 


12.5 ±0.50 


7.71 ± 1.00 


74- 


-76 


1.15 ±0.26 


1.07 ±0.26 


6-7 


10.4 ±0.50 


6.54 ±0.95 


76- 


-78 


1.69 ±0.25 


1.59 ±0.25 


7-8 


7.29 ±0.45 


4.79 ± 0.75 


78- 


-80 


0.73 ±0.22 


0.716 ±0.22 


8-9 


6.25 ±0.44 


4.13 ±0.73 


80- 


-81 


2.89 ±0.30 


2.88 ±0.30 


9-10 


5.37 ±0.40 


3.61 ±0.68 


81- 


-82 


2.05 ±0.28 


2.05 ±0.28 


10-12 


4.00 ±0.28 


2.68 ±0.45 


82- 


-83 


0.739 ±0.26 


0.739 ±0.26 


12-14 


3.93 ±0.30 


2.60 ±0.53 


83- 


-84 


1.42 ±0.26 


1.42 ±0.26 


14-16 


2.56 ±0.31 


1.66 ±0.45 


84- 


-85 


0.979 ±0.23 


0.979 ±0.23 


16-18 


2.86 ±0.35 


1.83 ±0.53 


85- 


-86 


1.3 ±0.23 


1.30 ±0.23 


18-20 


2.02 ±0.35 


1.28 ±0.47 


86- 


-87 


1.07 ±0.21 


1.07 ±0.21 


20-24 


2.66 ±0.28 


1.68 ±0.41 


87- 


-88 


0.861 ±0.21 


0.861 ±0.21 


24-28 


2.00 ±0.30 


1.25 ±0.39 


88- 


-89 


0.801 ±0.20 


0.801 ±0.20 


28-32 


1.98 ±0.31 


1.23 ±0.39 


89- 


-90 


0.444 ±0.18 


0.444 ±0.18 


32-36 


2.83 ±0.33 


1.74 ±0.46 


90- 


-91 


0.524 ±0.17 


0.524 ±0.17 


36-40 


2.61 ±0.33 


1.59 ±0.46 


91- 


-92 


0.479 ±0.17 


0.479 ±0.17 


40-44 


2.92 ±0.32 


1.78 ±0.49 


92- 


-93 


0.874 ±0.17 


0.874 ±0.17 


44-48 


2.85 ±0.31 


1.75 ±0.47 


93- 


-94 


0.345 ±0.14 


0.345 ±0.14 


48-52 


2.66 ±0.29 


1.69 ±0.43 


94- 


-95 


0.862 ±0.13 


0.862 ±0.13 


52-56 


2.72 ±0.28 


1.83 ±0.40 


95- 


-96 


0.380 ±0.11 


0.380 ±0.11 


56-60 


2.49 ±0.26 


1.77 ±0.35 


96- 


-97 


0.400 ±0.10 


0.400 ±0.10 


60-62 


2.56 ±0.35 


1.92 ±0.45 


97- 


-98 


-0.0245 ±0.080 


-0.0245 ±0.08 


62-64 


2.57 ±0.33 


1.99 ±0.43 


98- 


-99 


0.00468 ± 0.079 


0.00468 ± 0.07' 


64-66 


2.73 ±0.32 


2.18 ±0.41 


99- 


-100 


0.0227 ±0.13 


0.0227 ±0.13 


66-68 


2.80 ±0.316 


2.27 ±0.41 











Table A.6: Experimental double-differential cross sections at 20 degrees in the Ca(n,px) reaction for the neutron spectrum 
in table \4.3\ and for 94 MeV neutrons according to the subtraction method in section ed. 1.3\ 
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Bin 


Uncorrected 


Corrected 


Bin 


Uncorrected 


Corrected 


[MeV] 


[mb/sr/MeV] 


[MeV] 


[mb/sr/MeV] 


2.5-3 


7.64 ± 1.9 


4.59 ±2.2 


70- 


-71 


0.489 ±0.17 


0.428 ±0.19 


3-4 


11.3 ±0.57 


6.68 ± 1.0 


71- 


-72 


0.950 ±0.18 


0.834 ±0.22 


4-5 


6.43 ± 0.24 


3.87 ±0.51 


72- 


-73 


0.506 ±0.17 


0.456 ±0.18 


5-6 


12.6 ±0.34 


7.80 ±0.95 


73- 


-74 


0.459 ±0.16 


0.414 ±0.17 


6-7 


9.31 ±0.30 


5.84 ±0.78 


74- 


-75 


0.738 ±0.16 


0.689 ±0.17 


7-8 


6.97 ±0.28 


4.59 ±0.63 


75- 


-76 


0.355 ±0.15 


0.334 ±0.16 


8-9 


5.80 ±0.27 


3.84 ±0.60 


76- 


-77 


0.0924 ±0.14 


0.0861 ±0.14 


9-10 


5.16 ±0.25 


3.47 ±0.58 


77- 


-78 


0.700 ±0.15 


0.667 ±0.15 


10-12 


3.87 ±0.16 


2.60 ±0.38 


78- 


-79 


0.586 ±0.14 


0.571 ±0.14 


U-U 


2.72 ±0.15 


1.81 ±0.34 


79- 


-80 


0.366 ±0.13 


0.362 ±0.13 


14-16 


2.03 ±0.16 


1.33 ±0.31 


80- 


-81 


0.411 ±0.13 


0.409 ±0.13 


16-18 


2.29 ±0.17 


1.48 ±0.37 


81- 


-82 


0.0804 ±0.11 


0.0804 ±0.11 


18-20 


1.58 ±0.18 


1.01 ±0.31 


82- 


-83 


0.215 ±0.12 


0.215 ±0.12 


20-24 


1.47 ±0.14 


0.926 ±0.22 


83- 


-84 


0.274 ±0.11 


0.274 ±0.11 


24-28 


2.04 ±0.15 


1.28 ±0.29 


84- 


-85 


0.439 ±0.10 


0.439 ±0.10 


28-32 


1.78 ±0.15 


1.11 ±0.27 


85- 


-86 


0.172 ±0.094 


0.172 ±0.094 


32-36 


1.65 ±0.15 


1.02 ±0.27 


86- 


-87 


0.295 ±0.089 


0.295 ±0.089 


36-40 


1.48 ±0.15 


0.907 ±0.25 


87- 


-88 


0.0167 ±0.076 


0.0167 ±0.076 


40-44 


1.13 ±0.15 


0.694 ±0.22 


88- 


-89 


0.360 ±0.080 


0.36 ± 0.080 


44-48 


1.55 ±0.14 


0.971 ±0.26 


89- 


-90 


0.325 ±0.10 


0.325 ±0.10 


48-52 


1.22 ±0.13 


0.785 ±0.21 


90- 


-91 


-0.202 ±0.00 


-0.202 ±0.00 


52-54 


1.12 ±0.18 


0.746 ±0.27 


91- 


-92 


0.076 ± 0.074 


0.076 ±0.074 


54-56 


1.04 ±0.18 


0.720 ±0.25 


92- 


-93 


0.102 ±0.047 


0.102 ±0.047 


56-58 


1.26 ±0.17 


0.897 ±0.27 


93- 


-94 


0.0237 ±0.063 


0.0237 ±0.063 


58-60 


0.966 ±0.16 


0.712 ±0.23 


94- 


-95 


0.149 ±0.046 


0.149 ±0.046 


60-62 


0.851 ±0.16 


0.647 ±0.2 


95- 


-96 


0.0882 ±0.046 


0.0882 ±0.046 


62-64 


0.773 ±0.15 


0.605 ±0.18 


96- 


-97 


0.108 ±0.046 


0.108 ±0.046 


64-66 


0.871 ±0.14 


0.701 ±0.19 


97- 


-98 


0.0595 ± 0.044 


0.0595 ± 0.044 


66-68 


1.07 ±0.14 


0.876 ±0.19 


98- 


-99 


0.0190 ±0.045 


0.0190 ±0.045 


68-70 


1.07 ±0.13 


0.909 ±0.18 


99- 


-100 


0.0499 ±0.050 


0.0499 ±0.050 



Table A. 7: Experimental double-differential cross sections at 40 degrees in the Ca( n,px ) reaction for the neutron spectrum 
in table \4.3\ andfor 94 MeV neutrons according to the subtraction method in section \6.1.3\ 
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Bin 


Uncorrected 


Corrected 


Bin 


Uncorrected 


Corrected 


[MeV] 


[mb/sr/MeV] 


[MeV] 


[mb/sr/MeV] 


2.5-3 


0.0668 ±0.14 


0.0401 ±0.15 


30-35 


0.175 ±0.020 


0.114 ±0.082 


3-4 


5.41 ±0.42 


3.19 ±0.60 


35-40 


0.0721 ±0.015 


0.0472 ± 0.045 


4-5 


3.90 ±0.20 


2.35 ±0.35 


40-45 


0.0183 ±0.011 


0.0123 ±0.018 


5-6 


8.68 ±0.30 


5.37 ±0.69 


45-50 


0.00676 ± 0.0082 


0.00474 ±0.010 


6-7 


6.76 ±0.27 


4.26 ±0.59 


50-55 


-0.0278 ±0.0090 


-0.0208 ±0.031 


7-8 


4.21 ±0.21 


2.80 ±0.41 


55-60 


0.0146 ±0.0085 


0.0117 ±0.020 


8-9 


3.55 ±0.18 


2.37 ±0.39 


60-65 


0.0244 ±0.0099 


0.0209 ±0.033 


9-10 


3.37 ±0.16 


2.31 ±0.39 


65-70 


0.0274 ±0.0090 


0.0246 ±0.041 


10-12 


1.97 ±0.055 


1.36 ±0.20 


70-75 


0.00282 ±0.011 


0.00267 ±0.011 


12-14 


1.09 ±0.045 


0.762 ±0.15 


75-80 


0.0165 ±0.011 


0.0161 ±0.021 


14-16 


0.946 ± 0.046 


0.663 ±0.16 


80-85 


0.0562 ±0.013 


0.0561 ±0.014 


16-18 


0.709 ±0.052 


0.493 ±0.16 


85-90 


0.0346 ±0.015 


0.0346 ±0.015 


18-20 


0.283 ±0.053 


0.193 ±0.096 


90-95 


0.0914 ±0.016 


0.0914 ±0.016 


20-25 


0.420 ±0.024 


0.280 ±0.11 


95-100 


0.0274 ±0.019 


0.0274 ±0.019 


25-30 


0.149 ± 0.0233 


0.0968 ±0.057 









Table A.8: Experimental double-differential cross sections at 140 degrees in the Ca(n,px) reaction for the neutron spectrum 
in table \4. 3\ andfor 94 MeV neutrons according to the subtraction method in section \6.1.3\ 



Bin 


Uncorrected 


Corrected 


Bin 


Un-corrected 


Corrected 


[MeV] 


[mb/sr/MeV] 


[MeV] 


[mb/sr/MeV] 


2.5-3 


0.152 ±0.26 


0.0915 ±0.28 


30-35 


0.152 ±0.023 


0.0984 ±0.073 


3-4 


6.70 ± 0.60 


3.95 ±0.81 


35-40 


0.204 ±0.020 


0.133 ±0.12 


4-5 


6.44 ±0.38 


3.87 ±0.61 


40-45 


0.137 ±0.017 


0.0918 ±0.10 


5-6 


6.97 ±0.36 


4.31 ±0.62 


45-50 


0.0132 ±0.015 


0.00927 ±0.019 


6-7 


4.68 ±0.32 


2.95 ±0.50 


50-55 


0.0163 ±0.015 


0.0122 ±0.023 


7-8 


4.05 ±0.30 


2.69 ±0.46 


55-60 


-0.00418 ±0.017 


-0.00335 ±0.018 


8-9 


3.63 ±0.25 


2.43 ± 0.44 


60-65 


-0.00651 ±0.019 


-0.00558 ±0.021 


9-10 


2.89 ±0.21 


1.98 ±0.38 


65-70 


-0.0545 ±0.019 


-0.0491 ±0.081 


10-12 


1.59 ±0.083 


1.10±0.17 


70-75 


0.0539 ±0.021 


0.0511 ±0.073 


12-14 


0.915 ±0.079 


0.640 ±0.14 


75-80 


0.0449 ± 0.024 


0.0439 ±0.055 


14-16 


0.647 ±0.079 


0.454 ±0.14 


80-85 


0.134 ±0.023 


0.134 ±0.026 


16-18 


0.616 ±0.080 


0.428 ±0.16 


85-90 


0.146 ±0.028 


0.146 ±0.028 


18-20 


0.388 ±0.078 


0.266 ±0.13 


90-95 


0.254 ±0.030 


0.254 ±0.030 


20-25 


0.223 ±0.030 


0.148 ±0.063 


95-100 


0.229 ±0.036 


0.229 ±0.036 


25-30 


0.184 ± 0.0281 


0.120 ±0.070 









Table A.9: Experimental double-differential cross sections at 160 degrees in the Ca(n,px) reaction for the neutron spectrum 
in table \4.3\ andfor 94 MeV neutrons according to the subtraction method in section \6Jj\ 



48 



Bin 


Uncorrected 


Corrected 


Bin 


Uncorrected 


Corrected 


[MeV] 


[mb/MeV] 


[MeV] 


[mb/MeV] 


0-2.5 


_ 


_ 


46-48 


6.69 ±0.29 


4.15 ±0.41 


2.5-3 


12.1 ±0.32 


7.24 ±0.36 


48-50 


5.23 ±0.28 


3.32 ±0.37 


3-4 


102 ±0.52 


59.7 ±0.77 


50-52 


5.04 ±0.35 


3.26 ±0.55 


4-5 


73.7 ±0.26 


45.0 ± 0.47 


52-54 


4.54 ±0.33 


3.05 ±0.56 


5-6 


129 ±0.33 


78.7 ±0.73 


54-56 


4.46 ± 0.33 


3.07 ±0.49 


6-7 


97.3 ±0.30 


60.4 ±0.62 


56-58 


5.44 ±0.29 


3.77 ±0.53 


7-8 


69.7 ±0.26 


46.1 ±0.49 


58-60 


4.27 ±0.32 


3.12 ±0.44 


8-9 


59.2 ±0.24 


39.5 ± 0.47 


60-62 


3.85 ±0.31 


2.91 ±0.39 


9-10 


52.5 ±0.21 


35.4 ±0.45 


62-64 


3.68 ±0.30 


2.86 ±0.37 


10-12 


35.2 ±0.11 


23.6 ±0.26 


64-66 


4.00 ±0.27 


3.21 ±0.36 


12-14 


24.2 ±0.10 


16.4 ±0.23 


66-68 


4.54 ±0.29 


3.70 ±0.37 


14-16 


18.3 ±0.10 


12.0 ±0.22 


68-70 


4.50 ±0.28 


3.79 ±0.35 


16-18 


18.0 ±0.11 


12.0 ±0.26 


70-72 


3.25 ±0.24 


2.84 ±0.28 


18-20 


11.0 ±0.13 


7.19 ±0.22 


72-74 


2.53 ±0.27 


2.27 ±0.28 


20-22 


11.8±0.11 


7.34 ±0.24 


74-76 


2.34 ±0.45 


2.15 ±0.32 


22-24 


12.0 ±0.11 


7.35 ±0.27 


76-78 


2.23 ±0.25 


2.11 ±0.25 


24-26 


10.9 ±0.12 


6.82 ±0.28 


78-80 


2.42 ±0.12 


2.33 ±0.20 


26-28 


11.3 ±0.14 


6.95 ±0.30 


80-82 


2.89 ±0.56 


2.89 ±0.74 


28-30 


10.0 ±0.14 


6.25 ±0.31 


82-84 


1.55 ±0.23 


1.55 ±0.19 


30-32 


9.53 ±0.12 


5.99 ±0.31 


84-86 


1.71 ±0.15 


1.71 ±0.15 


32-34 


9.74 ±0.12 


6.04 ±0.39 


86-88 


1.29 ±0.21 


1.29 ±0.21 


34-36 


10.7 ±0.13 


6.64 ±0.46 


88-90 


2.47 ± 0.070 


2.47 ±0.070 


36-38 


8.17±0.11 


4.76 ±0.39 


90-92 


0.684 ±0.045 


0.684 ±0.041 


38-40 


9.51 ±0.12 


5.42 ±0.61 


92-94 


0.942 ± 0.044 


0.942 ± 0.044 


40-42 


5.76 ±0.13 


2.49 ±0.77 


94-96 


0.949 ± 0.060 


0.949 ±0.060 


42-44 


8.05 ±0.21 


3.59 ±0.68 


96-98 


0.516 ±0.15 


0.516 ±0.15 


44-46 


7.12 ±0.32 


3.94 ± 1.3 


98-100 


0.315 ±0.15 


0.315 ±0.15 


44-46 


7.12 ±0.31 


3.94 ±0.41 









Table A.IO: Experimental energy-differential cross sections in the Ca(n,px) reaction for the neutron spectrum in table \4.3\ 
and for 94 MeV neutrons according to the subtraction method in section ^. 1.3\ 



Angle 


Uncorrected 


Corrected 


[deg] 


[mb/sr] 


20 


253 ± 4.76 


175 ±6.69 


40 


162 ±2.51 


107 ±4.22 


140 


51.5 ±0.792 


33.7 ± 1.73 


160 


52.4 ± 1.14 


35.1 ±2.04 



Table A.ll: Experimental angular-differential cross sections in the Ca(n,px) reaction for the neutron spectrum in table \4.3\ 
and for 94 MeV neutrons according to the subtraction method in section ^. 1.3\ Cutoff energy is 2.5 MeV. 



Uncorrected Corrected 

^nb^ 

1220 ±3.21 788 ±5.83 



Table A. 12: Experimental total cross sections in the Ca(n,px) reaction for the neutron spectrum in table \4.3\ andfor 94 MeV 
neutrons according to the subtraction method in section \6.1.3\ Cutoff energy is 2.5 MeV. 
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B Notes about extracting information from NN databases 



When using databases to obtain data on elastic Nucleon-Nucleon (NNi scattering it is not always intuitive 



to get what one wants from it. The problem is that NN databases often gives the cross section for a given 
neutron scattering angle, but one most often knows the proton angle. To get the cross section one wants, some 
transformations between different frames and calculations of ratio-factors are needed [42] . This is an eight-step 
process to get through this without (too many) tears. First some necessary kinematic relations for a general 
scattering reaction ai^e given in 1.1 8J . 



mi 



Pi = dT( + 2Timi 



P3 



tan 6*3 



2^s 



TO2 + 2£'im2 
(s + TTig — 
Pi 

El + 1112 

{E^^Y - ml 



sm ( 



COSf 



(B.la) 
(B.lb) 
(B.lc) 
(B.ld) 

(B.le) 

(B.lf) 

(B.lg) 



The relevant quantities are defined in figure B.3 where is the mass, is the kinetic energy and is the 
momentum of particle i. Quantities in the CM system are labeled cm, unlabeled quantities are in the laboratory 
frame. 



P4 



Figure B.3: Scattering in the laboratory frame, adapted from [18% 



For a NN reaction with a neutron on a proton target we switch the indices in the notations as 1 , 4 ^ n and 
2, 3 ^ p. Known parameters are assumed to be m^, m„, T„ and 6*1^1^. An unlabeled 9 is always the scattering 
angle of the neutron. 

1 . Using the relations ( |B.la| i to ( |B.lg| i it is straightforward, although ( |B.lg| i might make it analytically tricky, 
to calculate 0™. 

2. The easy part: 6I™ = 180 - S™. 



3. Calculate 6'„ . This is preferably done via some computer code for relativistic kinematics, for example 
the code RELKIN circulating at Department of Neutron Physics, Uppsala University (INFi lf30l . 

4. From the |NN] database of choice, extract the cross section in the |CM| frame for the opposite case, that is 



when the neutron is scattered in the desired proton angle. That is, the cross section 



5. In the same way as in step|4] use a |NN| database to extract the cross section in the |CM frame for the 



desired case. That is, the cross section 
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Do stepj^l but in the laboratory frame, to get 



6 

7. Now to calculate the ratio factor, R, introduced in 



3n)I=e.a^fromthe 



NN 



database. 



for transformations between lab frame cross 



section and CM frame cross section. Since the R given in the reference only is valid for particles of equal 
masses the results from step|4jand step|6jis used to calculate R explicitly from the relation (3^)™ = 



,lab 



Using R obtained in step|7]it is now possible to transform the cross section from step|5]into the lab frame. 



by the relation previously used and obtained from |42|, (5^)g''_gii,b — (jS)™-^™ jj- 



Finally we have the cross section (g§)g'*giab which is just what we wanted, since the cross section of scat- 
tering a neutron at an angle of is equivalent to the cross section of scattering of a proton at an angle 0'*. 



Example Let us find the cross section for np scattering at 1 
the mass values r7i„ — 939.57 Me"' 
: 40.5 degrees and, consequentially. 



9[^[j — 20 degrees with an incoming neutron energy 
of T„ = 50 MeV. Using the mass values r7i„ = 939.57 MeV/c^ and to„ = 938.27 MeV/c^ the relations ( |B.la i 

qcm 



to (B.lg I will give 05, 



p 



dcr 

dcr 
dl7 



dCT 
df] 





w 69.4 degrees. 


cm 

1 




. 14.1"* 




=40.5 


sr 


cm 




mb 

=rf 15.0 — 


6'„ = 


139.5 


sr 


lab 




mb 

^ 53.9. — 


Bn = 


20.0 


sr 



139.5 degrees. Using RELKIN |30|, a code 
database |41 1 we get 



NN 



(B.2a) 



(B.2b) 



(B.2c) 

Note the different observable keywords DSG for cross section in the center of mass frame, and DSGL for the 



cross section in the laboratory frame. The ratio factor, R, can now be calculated using the information in ( B.2a 1, 
( B.2c| i and the definition in step|7] 



R 



dn/'e„=6i''''' 



14.1 
53i) 



0.261. 



(B.3) 



The ratio factor from ( |B.3[ ) and the cross section ( |B.2b| i can now be used to obtain 

daV" 

'l^/ e„ = 139.5 



daV 



which is the sought cross section 



0„=69.4 



1 = 15.0.^^57.5^^ 
R 0.261 sr 



(B.4) 
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